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We present an element selective resonant magnetic x-ray scattering study of NdFe3(BOs)4 as a 
function of temperature and applied magnetic field. Our measurements show that the magnetic order 
of the Nd sublattice is induced by the Fe spin order. When a magnetic field is applied parallel to the 
hexagonal basal plane, the helicoidal spin order is suppressed and a collinear ordering, where the 
moments are forced to align in a direction perpendicular to the applied magnetic field, is stabilized. 
This result excludes a non-collinear spin order as the origin of the magnetically induced electric 
polarization in this compound. Instead our data imply that magnetic frustration results in a phase 
competition, which is the origin of the magneto-electric response. 

PACS numbers: Valid PACS appear here 



INTRODUCTION 

The coupling between magnetism and electric polariza- 
tion in so-called multiferroic materials is a major topic 
of current condensed matter research, since it is of large 
interest for both basic science and technological appli- 
cations [1, 2]. Unfortunately ferroelectric and magnetic 
order rarely coexist and even if they do, the coupling 
between them is usually very weak [2]. 

It was therefore greeted with great excitement when 
frustrated magnetic materials were discovered, where the 
coupling of magnetic and ferroelectric orders is extraor- 
dinarily strong. Typically these systems relieve the mag- 
netic frustration to some extend via e.g. lattice distor- 
tions, which, as a byproduct, also create a ferroelectric 
polarization [1-3]. In this way the magnetic and ferro- 
electic (FE) orders are directly connected to one another, 
resulting in a very strong magneto-electric coupling. 

Typical examples are non-collinear magnets [4], where 
the spin order rotates about an axis e and propagates 
along a given direction q. In these magnetic phases a FE 
polarization P can exist, which, according to the conven- 
tional theory, is given by P oce x q [4, 5]. However, FE 
order can also be driven by collinear spin order [1, 6]. 
In order to reveal the microscopic origin of the ME cou- 
pling it is therefore very important to know the type of 
magnetic order. 

Among multiferroic systems that show non-collinear 



spin order, the magneto-electric (ME) NdFe3(BOs)4 has 
attracted considerable attention due to the large coupling 
between magnetic and electric orders [7, 8]. Upon cooling 
in zero magnetic field, a collinear magnetic order sets in 
first at Tn ~ 30 K, consisting of ferromagnetic hexagonal 
afr-planes, which are stacked antiferromagnetically along 
the perpendicular odirection (cf. Fig. 1) [9, 10]. Re- 
cent polarized neutron scattering experiments [10] have 
shown that a commensurate (C) to incommensurate (IC) 
transition is present at Tjc ~ 13.5 K. At this tempera- 
ture the magnetic structure turns into an IC spin helix 
that propagates along the c axis. Although this neutron 
study revealed important details about the helical spin 
order in NdFe3(B03)4, the microscopic origin of the spin 
helix remains to be clarified. 

Even more interesting, at low temperatures, the elec- 
tric polarization strongly increases when a magnetic field 
is applied parallel to a, reaching a maximum of P a ~ 
390 /iC/m 2 at B a ~ 1.3 T and 4.5 K[7, 11] (in TbMn0 3 
P c ~ 800 /iC/m 2 [3]). The spin helix in zero field cannot 
explain the electric polarization since e || q, for which 
the conventional theory yields P = 0. This already im- 
plies that the spin order in applied magnetic fields must 
change. As stated already above, it is essential to deter- 
mine these changes in order to uncover the driving force 
behind the ME coupling. Up to date, however, the spin 
order in applied magnetic fields was unknown for this 
compound. 
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FIG. 1. (Color online) (a) Crystal structure of NdFe3(BOs)4 with Fe chains of octahedra (brown), Nd prisms (gray) and the 
B triangles (green) . (b) Experimental scattering geometry showing the scattering plane (yellow) and the hexagonal basal plane 
(blue) where the a-axis is perpendicular to the scattering plane. Temperature dependent intensities of the magnetic (0,0,15/2) 
reflection at the Fe K (c) and Nd L2 (d) edges. The corresponding normalized integrated intensities are given in panels (e) and 
(f). Open circles in figure (e) show the T-dependent Fe moment, squared and scaled for comparison, from neutron diffraction 
(taken from [9]). Red scale in figures (c) and (d) show the (0,0,6) Bragg peak FWHM=0. 0018(3) r.l.u. Red lines are fits to the 
data as explained in the text. 



In this letter we clarify the two central questions men- 
tioned above, namely the origin of the spin helix as well as 
the magnetic order in applied magnetic fields. By means 
of resonant magnetic x-ray scattering (RMXS) we show 
that the magnetic order of the Nd-sublattice is induced 
by the Fe-spin order and determine the magnetic field 
induced changes in the spin order. The obtained data 
provide clear experimental evidence that the C-IC tran- 
sition and the ME effect are both related to the magnetic 
frustration between the Nd and the Fe sublattice. 

RMXS measurements at the Nd L 2 (6.726 keV) and Fe 
K (7.128 keV) edges were performed at beamlines ID20 
(ESRF in Grenoble, France [12]) and MagS (BESSY, 
Berlin [13]). A scheme of the scattering geometry used 
for the present experiments is illustrated in Fig. 1(b). The 
sample was mounted inside a cryomagnet with the hor- 
izontal scattering plane parallel to the 6*c plane of the 
sample. Measurements down to 2 K and with magnetic 
fields up to 2 T along the crystallographic a axis were per- 
formed, i.e., the magnetic field was aligned perpendicular 
to the scattering plane. During the measurements both 
the incoming and the outgoing photon polarizations were 
controlled, where a (a f ) and tt (tt') refer to linearly in- 
coming (outgoing) polarized light perpendicular and par- 
allel to the scattering plane, respectively. 



NdFe3(BOs)4 crystallizes in a hexagonal structure 
with space group R32 with the hexagonal unit cell pa- 
rameters a w 9.59 A, c w 7.61 A [9]. As shown in 
Fig. 1(a), the structure consists of chains of edge shared 
Fe06 octahedra along the crystallographic c axis. These 
chains are connected by BO3 triangles and NdC>6 prisms. 
NdFe3(BC>3)4 therefore hosts two coupled magnetic sub- 
lattices of, respectively, Fe and Nd sites. 

Below T^ ~ 30 K, the appearance of superlattice re- 
flections at q = (0, 0, 3/ ± |) (where I is an integer) re- 
veals a doubling of the unit cell along c (see Fig. 1(c) 
and 1(d)). This observation is in perfect agreement with 
recent neutron[9, 10] as well as x-ray scattering [14] re- 
ports. The simultaneous appearance of these reflections 
together with the magnetic order and the fact that the re- 
flected photons for incoming tt polarization are mostly a 
polarized, clearly identifies these reflections as magnetic 
superlattice peaks. 

Figs.lc and Id show high-resolution scans of the 
(0,0,15/2) magnetic reflection at different temperatures, 
which were measured at the Fe K and Nd L2 edge, re- 
spectively for B=0T. Below T^ ~ 30 K a single C- 
reflection (green circles) is observed, corresponding to 
the collinear AFM structure [9, 10, 14]. Upon further 
cooling down to Tjc ~ 16 K the peak splits into two 
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IC peaks at (0,0,15/2±5) due to the formation of a long 
period spin helix (black triangles) that propagates along 
the c direction[9, 10]. Surprisingly, although the C re- 
flection becomes very weak and broad, it never vanishes 
completely which indicates a coexistence of the C and IC 
phases at T < Tic- 

The value of S strongly increases below Tic as does 
the intensity of the IC reflections (open circles). This 
temperature dependent S reveals a continuous change 
of the helix periodicity. Specifically, we find that, in 
real space, the period of the spin helix decreases from 
-3900 A (-523 unit cells) at 14 K to -1123 A (-146 
unit cells) at 2 K. Again, these observations are in very 
good agreement with previous neutron results [10]. 

Within the kinematical theory of magnetic diffraction, 
the RMXS integrated intensity at the Fe K and the Nd L2 
edge is directly related to the ordered magnetic moment 
on the Fe and the Nd sublattice, respectively [15, 16]. In 
other words, RMXS enables to study the magnetic order 
on the Fe and Nd sublattices separately by tuning the 
photon energy to the corresponding absorption edge. 

Fig. 1(e) shows the temperature dependent integrated 
RMXS intensities measured at the Fe K edge and nor- 
malized to its value at the lowest temperature. As can be 
observed in this figure, the C peak first increases rapidly 
with cooling below Tat. Upon further cooling below Tic, 
the C reflection collapses whereas the IC reflections ap- 
pear and increase in intensity. 

A similar peak splitting behavior is also observed at the 
Nd L2 edge. However, here the temperature evolution of 
the intensities is very different. As shown in Fig. 1(f), 
the intensity enhancement of the C peak is much weaker 
in the range Tic < T <Tat. Below Tic the C peak 
intensity barely changes, while the IC peaks intensities 
increase pronouncedly. 

The blue circles in Figs, le and If correspond to the in- 
tegrated intensity of the (0,0,15/2) peak in the C phase 
plus the summed integrated intensities of the two satel- 
lites at (0,0, 1 5/2 ± S) . This total integrated intensity is 
proportional to (SFe) 2 at the Fe K and (SNd) 2 at the Nd 
L2 edge, where (SFe) and (SNd) denote the magnetic order 
parameter for the Fe and the Nd sublattice, respectively. 

The measured (SFe) 2 displays the typical behavior 
expected for a magnetic ordering transition and com- 
pares well to the corresponding neutron diffraction re- 
sults (open circles in Fig. 1(e)). 

(SNd) 2 exhibits a completely different temperature de- 
pendence and, in particular, does not show a clear phase 
transition at Tjy. To understand the magnetic ordering 
behavior of the Nd sublattice, the corresponding curve 
in Fig. 1(f) was fitted to a mean field model, which as- 
sumes that the Nd order is merely induced by the Fe 
ordering [17]: 
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FIG. 2. (Color online) Temperature dependence of the 
(0,0,15/2) peak in different applied magnetic fields. tt(j / and 
7V7v' scattering is shown in the left and right panels, respec- 
tively. The magnetic field was always applied parallel to the 
a axis. (g,h,i): Illustration of the spin directions as a function 
of magnetic field. Red arrows correspond to the easy axes of 
different magnetic domains in the c phase. 



In this equation ks is the Boltzmann constant, AFe-Nd 
is the effective Fe-Nd coupling constant, S^ d is the Nd 
ordered moment at T=0 K. Further, we used (SFe) 2 = 
(*Sp e (l — T/T/v)^) as a parameterization, where Sp e = 
5/2 is the ordered Fe moment at T = and /3 = 0.31(2) 
was determined by a fit to the experimental data (red 
line in Fig. 1(e)). The factor 3 accounts for the three 
Fe neighbors per Nd in the formula unit. Note that 
AFe-Nd (SFe) is the magnetic mean field due to the Fe 
neighbors of the Nd sites that induces the magnetic or- 
der on the Nd sublattice. 

The red solid line in Fig. 1(f) shows the fit of Eq. 1 to 
the experimental data. The best fit is obtained using 
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2.7 /i b determined from neutron diffraction [9] 
and AFe-Nd = 0.60(1) T//i#. The excellent agreement 
with the data implies that the magnetic order on the Nd 
sublattice is induced by the driving Fe spin order. 

The intensity maps shown in Fig. 2 display the temper- 
ature dependence of the magnetic (0,0,15/2) reflection 
measured at the Nd L2 edge in different applied mag- 
netic fields. According to the previous neutron scatter- 
ing studies, the magnetic moments of both Nd and Fe 



are always confined to the hexagonal ab plane [9, 10]. In 
this situation, the point symmetry of the Nd sites implies 
that the RMXS of the Nd sites can be described within 
the formalism given in Ref. 18 and 19, i.e., the intensities 
I n7r f and I^a' can be expressed as: oc (z a sin2#) 2 and 
lira' (— Zb* cos#) 2 , where is the scattering angle and 
z n are the projections of the magnetic moments along the 
crystallographic a and b* directions (cf. Fig. 1(b)). 

As can be seen in Figs. 2a and 2b, no qualitative differ- 
ence between the two intensities is observed for B=0 T 
apart from the fact that I n7r f / Ina' ~ 3. The relative in- 
tensities in the two scattering channels is the result of 
a particular domain structure. For instance, in the C 
phase the magnetic moments of the various domains are 
aligned along 3 equivalent directions within the hexag- 
onal ab plane (Fig. 2(g)). This yields z a ,z\>* ^ and 
intensities in both the hit' and the ircr' channel, where 
the value of 1^/1^' depends on the volume fractions 
occupied by the different domains in the volume probed 
by x-rays. The magnetic domains will not be discussed 
in the following. 

More importantly, there are already remarkable 
changes in the RMXS for 5=0.8 T (Fig.2(c) and 2(d)). 
First, the IC phase is suppressed and only observed at 
lower temperatures. Second, there is an overall inten- 
sity transfer from to ira' . This implies a reduc- 
tion of z a and an increase of 2:5*, i.e, the spins rotate 
as illustrated in Fig. 2(h). This observation agrees with 
magnetization[14, 20] and sound velocity[21] measure- 
ments which show that at 5=0.8 T the system is in the 
center of a hysteresis corresponding to a spin flop transi- 
tion. 

Increasing the magnetic field to B=2 T, the IC phase 
is completely suppressed, as can be seen in Fig. 2(e) 
and 2(f). In this field, pure ttct' scattering is observed, 
whereas the is very weak, amounting only to 0.5% 
of the 7rcr f intensity. This implies that z a = 0, i.e., the 
spin moments are now essentially perpendicular to B; 
i.e.; parallel to b* (Fig.2(i)). Note also that the intensity 
in Fig. 2(e) is equal to the sum of the intensities mea- 
sured in Fig. 2 (a) and 2(b), showing that at B=2 T a 
single collinear magnetic domain is formed. 

These field dependent measurements reveal that the 
maximum of magnetic field induced ferroelectric polar- 
ization P coexists with collinear spin order. The spin 
helix can hence be excluded as an origin for the magneto- 
electric effect in NdFe3 (1303)4. 

The microscopic insights obtained by RMXS enable to 
deduce a consistent scenario for the magneto-electric cou- 
pling in NdFe3 (1303)4, which is based on the magnetic 
frustration between the Fe and Nd sublattices: according 
to previous reports [9], Fe has a 3d 5 configuration, which 
implies that both the Fe-Fe and the Nd-Fe superexchange 
interactions are antiferromagnetic (AFM). Our data fur- 
ther shows that (Sp e ) 2 corresponds to the primary or- 
der parameter, whereas (SNd) 2 is merely induced. We 




(a) Fel 3d (b) 

FIG. 3. (Color online) (a) Crystal structure of 
NdFe3 (B03)4 showing the next nearest neighbor superex- 
change paths between Fe-Fe (blue) and Nd-Fe (green), (b) 
Magnetic frustration between Ndl and Fe2 ions. The 
strongest coupling is given by the intra chain Fe-Fe AFM cou- 
pling (blue). The Fel,3-Ndl (green) and Fe2-Ndl (black) are 
expected to be weaker. 

therefore consider the coupling of Nd to fully ordered Fe 
chains. As illustrated in Fig. 3, the AFM Nd-Fe interac- 
tions are frustrated, even though not all Nd-Fe couplings 
have the same strength. Upon decreasing temperature 
(SNd) 2 increases, i.e., an average magnetic moment on 
the Nd sites develops. This finite (SNd) 2 then also in- 
creases the energy related to the Nd-Fe frustration in the 
collinear phase. 

Without applied magnetic field, this Nd-Fe frustration 
can be relieved by forming a spin helix, which optimizes 
the magnetic interactions and gains exchange energy 
AEic by slight rotations Acf) of neighboring spins along 
c. Since (SNd) 2 increases upon cooling, AEic/A<p will 
also depend on temperature, which naturally explains the 
strongly temperature dependent period A ex 1/ Acj) of the 
spin helix. We therefore conclude that the spin helix is 
driven by the Nd-Fe sublattice frustration. 

By applying a magnetic field parallel to the afr-plane, 
the spin helix is rapidly suppressed and the spin rota- 
tions are no longer energetically favorable. We propose 
that within the magnetic field enforced C phase, the sys- 
tem choses an alternative way to relieve the magnetic 
frustration, which involves a modulation of the lattice 
structure and creates a FE polarization. We cannot de- 
termine these distortions based on the present data, but 
the 3-fold rotation axis of the space group R32 has to be 
broken in order to allow for a ferroelectric polarization 
parallel to the afr-plane. Interestingly, this removes the 
symmetry relations exactly between the Fe-sites shown 
in Fig. 3(a) [22]. The required symmetry reduction is 
therefore fully consistent with the deduced microscopic 
mechanism, which further supports our conclusions. 

NdFe3 (1303)4 hence belongs to a class of materials 
where the magneto-electric properties are driven by the 
magnetic frustration of two interacting subsystems and 
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where the degree of frustration can even be tuned by 
changing temperature. This frustration can be relieved 
either by forming a spin helix or by creating a FE po- 
larization. The FE collinear and the non-FE spin he- 
lix are therefore competing phases. Due to this compe- 
tition an external magnetic field can easily switch be- 
tween the two phases, causing the ME effect. We believe 
that a quantitative theoretical analysis of this mecha- 
nism may help to identify new magnetoelecric materials 
with optimized properties. The presented RMXS study 
also beautifully illustrates the power of this experimental 
technique and shows how this enables to unveil hidden 
microscopic physics that can hardly be accessed other- 
wise. This strongly motivates further studies of materi- 
als with different interacting magnetic sublattices using 
RMXS. 
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